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© Methode de determination des caracteristiques d'une formation souterraine produisant un f luide. 

@ L'invention a pour objet une methode de determination 

des caracteristiques physiques d'un systeme forme d'un puits 

et d'une formation souterraine contenant un fluids et com- 

muniquant avec Ie puits. On provoque un changement de 

debit dudit fiuide et on mesure une grandeur caracteristique 

de la pression P du fiuide a des intervalles de temps successlfe 

A t On compare en suite 

- d'une part, 1'evolution theorique du logarithme de la 
derivee P'o de la pression sans dimension en fonction du 
logarithme de to/Co, la derivee P'o etant par rapport a to/CD, tD 
representant le temps sans dimension et Cd I'effet de com- 
pression ou de decompression du fiuide dans le puits, avec 

- d'autre part, revolution experimentale du logarithme 
^ de la derivee A P' de la pression en fonction du logarithme des 

intervals de temps correspondents A t, la derivee A P' etant 
^ par rapport au temps t. On determine enfin de la comparaison 
(0 desdites evolutions theorique et experimentale le produit kh 
Y- de la permeabilite k par I'epaisseur de ladite formation h, et le 

coefficient C. 

If) 
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La presente invention concerne les essais de puits 
d'hydrocarbures, permettant de determiner les caracteristiques 
physiques du systeme forme d'un puits et d'une % formation 
souterraine (appelee aussi "reservoir") produisant des 
5 hydrocarbures a travers le puits. De fagon plus precise, 
1 T invention se rapporte a une methode selon laquelle le debit de 
fluide produit par le puits est modifie en fermant ou en ouvrant 
une vanne qui se trouve a la surface ou dans le puits. Les 
variations de pression resultantes sont mesurees et enregistrees 

10 au fond du puits en fonction du temps ecoule depuis le debut des 
essais, c f est-a-dire depuis la modification du debit. Les 
caracteristiques du systeme puits-f ormation souterraine peuvent 
etre deduites de ces donnees experimentales. Les donnees 
experimentales des essais de puits sont analysees en comparant la 

15 reponse de la formation souterraine a un changement de debit du 
fluide produit , avec le comportement de modeles theoriqiies ayant 
des caracteristiques bien definies et soumis au meme changement 
de debit que la formation etudiee. Habituellement , les variations 
de pression en fonction du temps caracterisent le comportement du 

20 systeme puits-f ormation et l 1 enlevement a debit constant de 
fluides, par l'ouverture d^he vanne dans le puits initialement 
ferme, est la condition d'essai qui est appliquee a la formation 
et au modele theorique, Lorsque leurs comport ements sont 
identiques, on suppose que le systeme etudie et le modele 

25 theorique sont identiques aussi bien au point de vue quantitatif 
que qualitatif. En d f autres termes, ces reservoirs sont supposes 
avoir les memes caracteristiques physiques. 

Les caracteristiques obtenues de cette comparaison dependent du 
modele theorique : plus le modele est complique, plus les 
30 caracteristiques qui peuvent etre determinees sont nombreuses. Le 
modele de base est represente par une formation homogene avec des 
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"""" """ rle " re « *>«rieure i-psr-ssbles st svec UD s 
extension -i 1 c • . - 

. . . * / £ * afiQxe - aeBifc *™ U formation est alors 
radial, dirige vers le puits. 

Cependant, le modele theorique le plus couramment utilise est 
Plus complique. II comp orte les caracteristiques du modele de 
base auquel on ajoute des conditions internes telles que 1'effet 

parietal ("skin effect" en anglais) et lUflW a 

. «*«giaxs; et 1 effet de compression ou 

ds depression du fluid, daus Is puits CVsll DO ts 6totag ... „ 
eusUis, L..„. t p. rIStal .„ dHlni pH _ Mef£lclmt 8 s ™ 
aractariss r.odo™.^ M „ ..^^ fc u « 
formation a d jacante ao puits . t .. (£K de C0Bpres6lM ou ae 
deco^ressiuu du fluids d„s Is puits sst caracterise p ar _ 

produit psr 1. puits. eutre Is formation soutsrrsins st Is tets 
de puits. lots,u. use v.nn. sltuae en tSte de ^ 
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veria ions ds ptession. „ Wtlon du £emps> du ^ 

du Puits. ces coursas scut habitually ttscees s„ ooordonnsss 

dil " " " 5ChCUe 1 °«' rlttaI '°«. «• 

etsut port, su sbseisss. Ds plus. chs,„s course est cstsctstises 

cluT PlUSleUIS nratt " ~ ^sentent 
chs™ ™ «««etisti q us (ou uos cosMnsiso. de 
car.ctsristi,uss).du systsae thsori,ue for.a pst uu puits st „„ 
ressrvoir. On psts»itre s.ns dimension set dMini psr le 
psts„etts tie! (l a p« sslon p „ .„„ pl . ) 

expression ,ui indut cettsiues csract«ri,ti,„s, du systtae 
puits-rsssrvoir ds fs„n . randra is p aranetra sa „ s dime „ SIO „ 
indepsndant ds css caracteristiouss. 
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C'est ainsi que le coefficient S caracterise uniquement I'effet 
parietal mais est independant des autres caracteristiques da 
reservoir et des conditions experimentales telles que le debit, 
la viscosite du fluide, la permeabilite de la formation, etc. 

5 Lorsque le modele theorique et le systeme etudie puits- format ion 
correspondent, la courbe experimental et l'une des courbe-types 
representees avec les memes echelles de coordonnees ont une forme 
identique mais sont decalees I'une par rapport a I'autre. Les 
decalages suivant les deux axes, en ordonnee pour la pression et 

10 en abscisse pour le temps, sont proportionnels a des valeurs de 
caracteristiques du systeme puits-reservoir qui peuvent ainsi 
etre detenainees. 

Des informations qualitatives sur la formation souterraine, telle 
que par exemple la presence d'une fracture, sont obtenues par 

15 l f identification des differents regimes sur le graphs en echelle 
logarithmique representant les donnees experimentales. Sachant 
qu'une caracteristique particuliere du systeme puits-reservoir, 
telle que par exemple une fracture verticale, se caracterise par 
un regime particulier, tous les differents regimes apparaissant 

20 sur le graphe des donnees experimentales sont identifiers pour 
selectionner le modele de systeme puits-reservoir approprie. Des 
graphes specialises ne tenant compte que d'une partie des donnees 
experimentales permettent la determination plus precise des 
caracteristiques du systeme. On revient ensuite an graphe en 

25 echelle logarithmique prenant en compte toutes les donnees pour 
confirmer le choix du systeme et la determination quantitative 
des caracteristiques de la formation. Ces dernieres sont 
obtenues en selectionnant une courbe-type ayant la meme forme que 
la courbe experimental et en determinant le decalage des axes de 

30 coordonnees de la courbe experimental par rapport a la courbe 
theorique. ... 



0125164 



A un meme modele theorique correspondent plusieurs graphes de 
courbes- type. Cecf. depend de, pa Tare s r . e? sans d««-.«- 
pour la representation des axes de coordonnees du graphe. ain^i 
que d'un ou plusieurs "index". On "index" n'est autre qu'un 
parametre supplfetaire (ou une co m binaison de parametres)choisi 
pour la representation des courbes, en co.ple.ent des parages 
sans dimension des axes de coordonnees. La comparison des 
differentes methodes utilisees est donnee dans 1'artide intitule 
A Comparison Between Different Skin and Wellbore Storage Type 
10 Curves for Early-Time Transient Analysis" (Comparison entre les 
differentes courbes-type avec effet parietal et effet de 
compression ou de decompression pour 1'analyse des transitoires 
precoces) par A.C. Gringarten & al., publie par "Society of 
Petroleum Engineers of AIKE", (n°SPE8205). Le brevet des 
ISEtats-Onis d'Amerique 4,328,705 decrit egalement une methode 
selon laquelle les courbes-type sont representees en utilisant la 
pression sans dimension P D pour 1'axe des ordonnees et le rapport 
V C D I'axe des abscisses, t, etant le temps sans dimension 

et C D le coefricient sans dimension caracterisant 1'effet de 
20 compression ou de decompression du fluide dans le puits 
L inconvenient de la methode decrite dans ce brevet est que les 
courbes-type ont des formes variant relativement lentement 1W 
Par rapport a 1'autre. U en resulte une certaine incertitude 
dans le choix de la courbe-type correspondant a la courbe 
25 experimental On remar q Ue egalement que pour effectuer une 
analyse complete, on est oblige de faire appel non seulement a un 
graphe en echelle logarithmique representant 1' ensemble des 
donnees experimental , mais d'utiliser egalement des graphes 
specialises, en echelle semi-logarithmioue par exemple, pour 
30 n' analyser gu 1 ™ _ , 

^ — P^.xe aes donnees mais d'une fagon plus 

precise. 
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On a dija songe a utiliser la derivee raathematique de la presslon 
sans dimension, P' D , au lieu de la pression sans dimension P D< 
C'est ainsi que dans 1' article ' intitulg "Application of the P'^- 
Function to Interference Analysis" (Application de la fonction 
5p' D i l'analyse des tests d' interference) publie dans "Journal of 
Petroleum Technology" d'aout 1980, page 1465, Involution de la 
derivee P' D (derivee par rapport a tp) en fonction de est 
utilisee pour analyser des essais d' interference entre un puits 
de production et un puits d' observation. On enregistre les 

lOvariations de pression dans le puits d' observation lorsqu'on 
modifie le debit de fluide produit par le puits de production. 
Dans ce cas, l'effet parietal et l'effet de compression ou de 
decompression du fluide n'interviennent pas. C'est done un cas 
tres simple dans lequel on analyse la reponse de la formation 

I5 S outerraine dans un puits eloigne du puits producteur. II en 
resulte qu'il n'y a pas une famille de courbes-type mais une 
seule courbe. 

La derivee de la pression P' D (derivee par rapport a t p ) a ete 
egalement mise a profit pour caracteriser les reservoirs 

20 C ontenant deux failles etanches paralleles autour du reservoir, 
dans 1' article intitule "Detection and Location of Two Parallel 
Sealing Faults Around a Well" (Detection et localisation de deux 
failles etanches paralleles autour d'un puits) paru dans "Journal 
of Petroleum Technology" d'octobre 1980, page 1701. Cet article 

25ne traite que d'un probleme particulier. 

Le comportement en pression d'un puits produisant un fluide 
legerement compressible a. travers un plan unique d'une fracture 
verticale dans un reservoir infini a ete analyse a l'aide de la 
derivee mathematique de la pression sans dimension P' D> (derivee 

30par rapport a un temps sans dimension t_,) dans 1' article 

- ' • ■ - Df 

• intitule "Application of P' D " Function "to Vertically Fractured 
Wells" (Application de la fonction P' D au puits fracture 
verticalement) paru dans "Society of Petroleum Engineers" of 
AIME, SPE 11028, du 26-29 septembre 1982. 
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Cet article ne coucerne qu'un cas particulier dans l equ el la 
couroe-ty pe e6fc UIllque t i.„„„i , 

, P Ur le <l uel les avantages de 

1 atilisatioa de la dirivaa da U prasslan ne saar pas aLaats 
par rapport ,»* „ 6thooes classl , ues . ^ ^ 

laffet de saapr.ssion au da decompressiaa d„ fl„ lde 
n interviennent pas. 

U presents iaventien a pa„r aajec aaa -echede da detertrtnatiea 
des earae.eristieeas d .„„ systSme pults . r6se „ oll 

"-ificatias aatra le ce.p„rte„ent eKperi.ental d„ 

at la ceaperteaeat d-ua ™ dHe th«arl,aa. c. aadtte .st gaaatal 
a sav.ir ,ae 1. f OImLtioa DMt ^ fc ^ 

aeil tiaat eeapte da Verfet parietal at da l.. ££et de 
eeapressian „» da deceapresaie. da fl ul da at iva.taallaa.at da la 
"daabla parasite du reservoir at das fractaras da p„i ts . U 
aethoda seloa la prasaata iavaatlaa peraet aaa analys. 8 l ot ala at 
aaroue da ecnperteaeat da svsteae pnits-reservair, saas raaourir 
a_ das aaalvses speeiaUsees. Wnveatien peraet eealeaent 
1 eealvse das deanees experi-eatales lars,aa 1, caaditiaa iapasee 
an systaaa ast la feraetnre du pairs, grSce , _ choix ^ 
das paraaatraa. la airbeds salaa la prasaata inveatioa paat 
e 8 ale«*„t atra avaatageaseaent eeabinie i „ne aethede da 1'art 
anterieur. 

De facon plus precise, la presente invention concerne une methode 
5de detennination des caracteristiques P h ys i ques d'un system 
forme d un puits et d'une formation souterraine ccntenant un 
fluxde et communiquant avec ledit puits, ladite formation 
presentant de 1'effet parietal et/ou le fluide se co.pri.ant ou 
se deco.pri.nant dans le puits et ladite formation etant homogene 
: hgt5 ^ Sne - Selon la methode, on provoque un changes de 
debit dudit fluide, on mesure une grandeur caracteristique de la 
Pression P du fluide a des intervalles de temps successifs At et 
on compare, 
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d'une part, a partir d'un modele theorique de systeme de 
puits-reservoir , 1' evolution theorique du logarithme de la 
derivee P'^ de la- pression sans dimension en fonction du 
logarithme de t D /C D , ladite derivee etant par rapport a 
5 S'S* *D re P r5sentant le temps sans dimension et C D le 

coefficient sans dimension de I'effet de compression ou de 
decompression du fluide dans le puits, avec 
d f autre part, Involution experimentale du logarithme de la 
derivee i P' de la pression en fonction du logarithme des 
10 intervalles de temps correspondants At, ladite derivee A p 1 

etant par rapport au temps t, 
et on determine de la comparaison desdites evolutions theorique 
et experimentale au moins une caracteristique du systeme 
puits-formation, choisie parmi le produit kh de la permeabilite k 
15 par I'epaisseur de ladite formation h, le coefficient et le 
coefficient S de I'effet parietal, 

Ladite evolution theorique peut avantageusement etre celle du 
logarithme du produit P' .t D /C D en fonction du logarithme de 
t^/Cjj et ladite evolution experimentale est celle du logarithme 
20 du produit A P f . At en fonction du logarithme de At. 

Ladite evolution theorique peut egalement etre fonction d'un 

2S 

index representant une grandeur caracteristique du produit C^e . 
Lorsque le changement de debit du fluide correspond a la 
fermeture du puits, on peut avantageusement comparer ladite 
25 evolution theorique avec devolution experimentale du logarithme 
de l f expression : 

t + A t 

.At . A P • 

t 

P 

en fonction du logarithme des «intervalles. .de-temps. A, t> 9 +t^ etant 

p 

le temps pendant lequel le puits a ete mis en production. 
30 Certaines e tapes de la presente invention, notamment 
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1' identification des donnees experimentales avec le comportement 
d ; un modele theorique ayant des caracteristiques bien precises, 
peuvent etre mises en oeuvre a l'aide d'un ordinateur. Cependant, 
ces Stapes sont avantageusement mises en oeuvre en tracant un 
5 graphe theorique en coordonnees cartesiennes et en echelle 
logarithmique, ledit graphe representant Involution theorique de 
la derivee P« D en fonction de ou encore revolution 

theorique du produit P'^VS en fonction de t^Cp. 
On peut egalement tracer une courbe experimental a l'aide des 

10 donnees experimentales avec la meme echelle logarithmique que 
ledit graphe theorique, la courbe experimental representant 
soit Involution experimental de A P* en fonction de At, soit 
Involution experimental du produit A P'. At en fonction de 
A t. On peut alors faire correspondre la courbe experimental 

15 avec l'une des courbes-type du graphe thgorique et on peut en 
determiner certaines caracteristiques physiques du systeme 
puits-formation souterraine. 

L« invention a egalement pour objet les graphes theoriques obtenus 
tel qu'indique precedemment. 
20 L' invention sera mieux comprise a l'aide de la description qui va 
suivre de modes de realisation de l'invention donnes a titre 
d'exemples explicatifs mais nullement limitatifs. La description 
se rapporte aux dessins qui l'accompagnent dans lesquels : 

la figure 1 represente en echelle logarithmique un graphe de 
25 courbe-type representant P' D en fonction de tp/C^ l'index 

representant les valeurs de C D e 2S ; 

la figure 2 montre un graphe de courbes-type en echelle 
logarithmique representant P'.t^ en fonction de t_/C , 
l'index etant C D e 2S ; * D 

30 - la figure 3 illustre la methode selon la presente 

invention pour la determination des caracteristiques 
physiques d'une formation souterraine produisant un 
fluide ; 
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la figure k represente en echelle logarithraique un graphe de 
courbes-type representant P'p.tp/Cp en fonction de t^/C^ 
pour une formation souterraine a double- porosite ; et 
la figure 5 represente deux series de courbes-type en 
5 echelle logarithmique , l'une raontrant des courbes-type 

connues et 1' autre montrant des courbes-types selon la 
present e invention, 
Avant de mettre un puits d 1 hydrocarbures en production, on 
effectue generalement des mesures de fagon a determiner les 
10 caracteristiques physiques de la formation souterraine produisant 
ces hydrocarbures. Cette etape preliminaire avant la production 
est tres importante car elle permet de definir les conditions les 
plus appropriees pour produire ces hydrocarbures et 
eventuellement pour ameliorer la production. L'une de ces mesures 
15 consiste a faire varier le debit du fluide produit, en ouvrant ou 
f ermant une vanne placee en tete de puits ou dans le puits 
lui-meme, et a enregistrer les variations resultantes de pression 
en fonction du temps ecoule a partir de la modification de debit 
du fluide produit. On peut par exemple fermer completement le 
20 puits et enregistrer 1' augmentation de pression resultante (on 
obtiendra alors une courbe experimentale dite courbe de 
"build-up"). On peut aussi faire produire un puits dont la 
production avait ete arret ee et enregistrer la diminution de 
pression correspondante (la courbe experimentale obtenue 
25 s'appelle alors "courbe de drawdown"). 

Les variations de pression en fonction du temps peuvent etre 
suivies a l'aide d'une sonde descendue dans le puits a 
l f extremite d'un cable. Le cable peut etre electrique et dans ce 
cas transmettre les informations de pression directetnent a un 
30 enregistreur situe a la surface. Lorsque le cable n'est pas 
conducteur, les variations de pression sont enregistrees dans des 
memoires disposees dans la sonde. Ces memoires sont ensuite lues 
en surface. On peut egalement installer une jauge de pression 



10 
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dans une poche laterale de la colonne de production du puits, a 
pr ° SiBitS de U fcras£loa Patrice. Uu .able conducts situe 
dans l'espace annulaire compris entre la colonne de production et 
le tubage relie la jauge de pression jusqu'a un enregistreur 
situe a la surface. Un tel dispositif est decrit par exemple dans 
les brevets des Etats-Unis n° 3,939,705 et 4,105,279. 
Les valeurs mesurees par les sondes de pression ne correspondent 
gineralement pas a la press ion elle-meme, mais S une grandeur 
caracteristique de la pression telle que , par exemple, une 
•0 difference de deux frequences. Par la suite, on utilisera 
l'expression "valeur de la pression" par coamodite et clarte 
tout en gardant a 1'esprit que les donnees experimental peuvent 
correspondre a une grandeur caracteristique de la pression. 
La figure 1 represente un graphe de nouvelles courbes-type, en 
15 gchelle logarithmique, representant la derivee mathematique P' 
de la pression sans dimension P D en fonction du rapport t /C D 
t D represente le temps sans dimension et C]) representant la 
valeur du coefficient sans dimension de 1'effet de compression ou 
de decompression du fluide dans le puits. La derivee mathematique 
20P' D est faite par rapport a t^. De plus, les variations de la 
derivee de^la pression P' D SO nt representees par rapport a un 
index C D e , qui n'est autre qu'une combinaison de deux 
caracteristiques physiques Cj) et S du systeme puits-reservoir 
analyse. 11 est a remarquer que 1'index C D e 2S peut prendre une 
25 valeur quelconque, qui n'est pas necessairement entiere. La 
valeur de la pression sans dimension P D est donnee par V equation 
suivante, en employant le systeme d'unites utilise couramment 
dans 1« Industrie petroliere et appele "Oil field units" page 185 
du livre intitule "Advances in Well Test Analysis" publie par 
30 "Society of Petroleum Engineers of Aih£ !: - 1977 : 
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dans laquelle : 

k represente la permeabilite de la formation souterraine, 
h'esr l'gpaisseur *de* la formation ; . / . . 

A P est la variation de pression, 
5q est le debit du fluide en surface, 
B est le coefficient de dilatation du fluide entre le reservoir 
et la surface (appele en anglais "formation volume factor") et 
is est la viscosite du fluide. 

La dgrivee mathematique P 1 ^ de la pression sans dimension P^, par 
10 rapport a t^/C^, est donnee par 1' equation suivante : 

p 1 = C ^pi 

D 0,04194 qB (2) 

dans laquelle A p f est la derivee (par rapport au temps t) de la 
variation de pression A p en fonction de I'intervalle de temps 
At qui represente l ! intervalle de temps s'etant ecoule depuis le 
15 debut de 1'essai de la formation, c'est-a-dire I'intervalle de 
temps entre V instant de mesure et 1' instant de modification de 
debit du fluide. 

La valeur du rapport t D /C D > dans le meme systeme d 1 unites que 
pour les equations precedentes, est donnee par : 

20 0,000295 kh . At 

C D ™ y C (3) 

dans laquelle C est l'effet de compression ou de decompression du 
fluide dans le puits. 

Le graphe de la figure 1 caracterise le comportement d f un modele 
de reservoir homogene et d'un puits presentant l'effet parietal 
25 et l'effet de compression ou de decompression du fluide dans le 
puits. 

Ce graphe £6t^Qbteriue>a> partir -de £ I 1 equation (A. 2) .de f'airticle 
intitule "Determination of fissure volume and block size in 
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fractured reservoirs by type-curve analysis" (Determination du 
voxume des fissures et de ia taiiie des blocs dans les reservoirs 
fissures par analyse par courbes-type) publie par "Society of 
Petroleum Engineers" en septembre 1980, n» SPE 9293. Cette 
equation est donnee dans le domaine de Laplace. Inversion dans 
le domaine des temps reels est obtenue a 1'aide d'un algorithme 
d inversion, tel que celui decrit par exemple par H. Stebfest 
dans Communications of the ACM, D-5" du 13.01.70, n*l page 47. 

Les courbes de la figure 1 se caracterisent par trois parties 
dxstinctes : la .partie gauche du graphe correspondent an* temps 
courts et etant caracteristique de 1'effet de decompression du 
fluide du puits (cet effet est le plus important a 1'ouverture de 
la vanne) ; ia partie droite du graphe corresponQant . ^ 

ecoulement pur radial du reservoir et une partie intermediaire 
entre les parties gauche et droite correspondant I un regime 
d ecoulement transitoire entre les deux ecoulements-limite 
precedents. Cet ecoulement intermediaire est fonction de 1'effet 
de decompression du fluide et de 1'effet parietal. 
Sur la partie gauche du graphe, les courbes tendent vers une 
asymptote correspondant a une derivee ggale a 1. En effet, au 
tout debut des essais, le phenomene predominant est 1'effet de 
decompression du puits lequel est caracterise par 1'equation : 
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La derivee de la pression sans dimension par rapport a t^/C^ peut 
s'ecrire : 



g (Pp) 

d <W 



P ' D - 1 (5) 



On voit que la derivee P» D , pour ce type d'ecoulement est egal a 
5 1 et que les courbes-type se reduisent a une droite de pente 
nulle. La partie droite du graphe de la figure 1, qui correspond 
a un ecoulement radial infini dans une formation homogene, est 
caracterisee par 1' equation : 

P D = 0,5 £ln Id + 0,80907 + In C D e 2S j 

10 In representant le logarithme neperien. 

En derivant Pp par rapport a t^/C^, on obtient : 



(6) 



d < VV " t D /c D 

et en passant en echelle logarithmique : 
15 log P f D = log 0, 5 - log ±JL 



(8) 



On remarque que la courbe representee par 1' equation (8) est une 
droite de pente egale a -1. Pour les temps courts et les temps 
longs, les courbes sont rectilignes et independantes de C D e 2S , ce 
qui est un avantage considerable par rapport aux methodes de 
20 l f art anterieur. Entre les deux asymptotes, pour les temps 
intermediates, chaque courbe d' index C n e 2S a une forme 
differente bien contrastee. 

Si dP represents la difference de deux mesures successives de la 
pression du fluide dans le puits et si dt represente l f intervalle 
25 de temps (court) separant ces deux mesures successives, on 
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calcule pour tous les couples de mesures successives les valeurs 
* r = dP/dt. Ce caicul permet de determiner de fa C on pratique 
les valeurs successives de la derivee mathematique AP' qui par 
definition est egale au rapport dP/dt lorsque dt tend vers 
Szero.En tragant la courbe A P' en fonction de At (At etant 
l'intervalle de temps entre l'instant de la mesure consideree et 
l'instant de modification de debit du fluide) de fagon a former 
un graphe experimental, et en prenant les memes echelles 
logarithmiques que celles adoptees pour tracer les courbes-type 
10 de la figure 1, on peut determiner les caracteristiques physiques 
du systeme puits-formation souterraine. En effet, le decalage des 
ordonnees de la courbe experimental et des courbes-type permet 
de determiner la valeur de C (ce qui est evident au vu de 
l'equation (2), en effectuant log P' D -logAP' et en connaissant 
les valeurs de q et B). Le decalage des abscisses de la courbe 
experimental par rapport a la courbe-type choisie permet de 
determiner la valeur kh (connaissant C et p, ce qui est evident 
au vu de I' equation (3) en effectuant logt D /C D - log At). Enfin 
le choix de la courbe-type correspondent a la courbe 
20 experimental permet de determiner le coefficient S (par le 
caicul prealable de C]) a partir de 1« equation (14) comme il sera 
montri par la suite) . Le graphe theorique de la figure i etant 
utilise de la meme fagon que celui de la figure 2, par 
comparison avec la courbe experimentale, seule l'utilisation du 
25 graphe de la figure 2 est illustree (figure 3). 

La methode de determination des caracteristiques physiques par 
utilisation du graphe de la figure 1 a ete amelioree en suivant 
devolution, non plus de la derivee mathematique de la pression 
^ sans dimension, mais en . suivant Involution, en fonction de 
J ° W du P roduit de ^ derivee P' D de la pression sans dimension 
(derivee par rapport a t^) par le rapport t^C.. Cette 
nouvelle methode est illustree a la figure 2 par un graphe 
representant le comportement d'une formation homogene presentant 
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l'effet parietal et l'effet de decompression du fluide du puits. 
L'axe des ordonnees correspond a P^.tp/C et l'axe des abscisses 
correspond a t D /C D> Pl D etant la derivee.de P D par rapport a 

2 S 

5 De plus, 1 1 index C D e a ete choisi pour representer les 
courbes-type. Comme dans le cas de la figure 1, l'effet 
predominant au debut de I'essai de puits est l'effet de 
decompression du fluide. Cet effet correspond aux equations (A) 
et (5). A partir de l 1 equation (5), on peut ecrire : 

D W CD ' } 

On remarque dans cette derniere equation, que pour les temps 
courts, les courbes-type tendent vers une asymptote de pente 
egale a 1. 

Pour les temps longs, correspondant a la partie droite du graphe 
15 de la figure 2, les equations (6) et '(7) demeurent valables 
puisqu'on se trouve en fin de test en ecoulement radial infini 
pour une formation homogene. L' equation (7) peut s 1 ecrire : 

P 'd • ^ °' 5 (10) 



D 



II en resulte que pour les temps longs, la valeur du produit 

20 pl D #t D^ C D eSt g 8 ale S °» 5 et ^ ue les courbes-type tendent vers 
une asymptote de pente nulle. 

On remarque que pour le regime d' ecoulement intermediaire situe 
au centre du graphe de la figure 2, les courbes-type sont de 
forme tres contrastee, ce qui permet une identification beaucoup 
25 plus precise de la courbe experimentale avec I'une des 
courbes-type que par les methodes de l'art anterieur. Par rapport 
au graphe de la f igure. 1 , ,on peut dire que le . graphe. de la .figure 
2 correspond, en premiere approximation, k une rotation de 45° du 
graphe de la figure 1. Cependant, les courbes-type ont un relief 



0125164 



16 



Plus accents et l a presentation du graphe de la figure 2 est 
Plus pratique.Les valeurs de r index % f* sonc lndiquges sur ^ 
courbes-type. La figure 3 illustre 1'utilisation du graphe des 
courbes-type de la figure 2. Ce g rap he a ete reproduit sur la 
figure 3 avec, en ordonnee, F VVS et en abscissa t/n. Les 
differences de pression dP usuries dans le puits, pour des 
differences de temps successifs dt, sent utilises pour calculer 
les valeurs A P - « dP/dt conme indiqug prgcgdenment> ^ 
les valeurs successives de AP' par les int ervalle S de temps 
correspondents A t et on trace ensuite un graphe experimental 
representant le produit AP '. At en ordonnge eQ de & t 

en abscisse. Les valeurs de A P sont en psi (1 psi . ^ bar) 
et les valeurs de At sont en heure. Les deux g raphe s thiorique 
et experimental ont la meme ichelle logarithmique. On commence 
par superposer la partie droite, qui est rectiligne, de la courbe 
experimental tracee sur la figure 3 a 1'aide de points, a l a 
partie rectiligne des courbes-type a droite du graphe. Ceci est 
facile a realiser puisque cette partie des courbes est une droite 
de pente nulle. On decale ensuite le graphe experimental le long 
de 1 axe des temps de f aC on a faire corresponds sa partie gauche 
avec la partie gauche des courbes-type. Ceci est egalement facile 
puisque cette partie des courbes-type est une droite de pente 
egale a 1. Si la formation souterraine etudiee a un comportement 
hoogene, la courbe experimental doit alors se superposer 
parfaitement, aux erreurs de precision des mesures pres. a une 
courbe-type. Dans 1'exemple montre sur la figure 3, cette 
courbe-type correspond a C D e 2S = 10 10 . Le decalage des axes de 
coordonnees de la courbe experimental avec les axes des 
» courbes-type permet de determiner les valeurs du produit kh et la 
3^ valeur de 1'effet de depression du fluide dans le puits. En 
effet, en combinant les equations (2) et (3),on obtient : 



0125164 



17 



° C Q 141,2 qBy (ID 

qui s'ecrit : 

log (P' _M- log(AP'.At) = log kh (12) 

\ V 141,2qBy 

Le membre de gauche de cette derniere equation correspond au 
5 decalage des ordonnees represente par Y sur la figure 3. 
La valeur de Y pennet de determiner le produit kh. En effet, la 
valeur du debit de fluide q est generalement connue par des 
mesures ayant ete prealablement effectuees avec un debitmetre ou 
un separateur, et les valeurs du coefficient B de dilatation du 
lOfluide et de sa viscosite u sont determinees par 1' analyse 
d'echantillons du fluide (analyse appelee habituellement "PVT"). 
En consequence, la valeur du produit permeabilite - epaisseur kh 
peut etre determinee en connaissant la valeur Y mesuree. 
De la meme facon, 1* equation (3) peut s'ecrire 

15 log la - log At = log 0,000295 ^ (13) 
C D ^ c 
Le membre de gauche de cette equation correspond au decalage X 
des abscisses de la courbe-type choisi et de la courbe 
experimentale. Connaissant la valeur de ce decalage X ainsi que 
les valeurs de la viscosite p et du produit kh, on deduit de 

201' equation (13) la valeur du coefficient C de l'effet de 
decompression du fluide du puits. 

La valeur du coefficient S de 1' effet parietal est determine par 

la correspondance de la courbe experimentale avec l'une des 

courbes-type, la correspondance des deux courbe* conduisant a la 
2S 

v 25yaleur de C D e . La. .valeur de C D est determinee .par -la valeur 'de 

C par 1' equation suivante : 
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r = 0.8926 C 

0 °t a (1A) 



dans laquelle 0c t h represente le produit porosite-compressibilite- 
epaisseur, connu a partir d'etudes geologiques (telles que par 
5 exemple, l'analyse d'echantillons ou les diagraphies electriques) 
et r est le rayon du puits. La valeur du coefficient S pent done 
etre calculee a partir de la valeur de C D e 2S . 

Les courbes-type representees sur les figures 1 et 2 
correspondent au comportement d'un modele theorique de formation 

10 homogene lorsqu'on impose soudainement une augmentation du debit 
du fluide produit par la formation, et plus specialement 
lorsqu'on ouvre une vanne en surface du puits pour le faire 
produire a debit constant alors qu'il etait ferme auparavant 
(courbe de "drawdown"). 

15 Selon 1'une des caracteristiques de la presente invention, pour 
les analyses d'essais de puits correspondant a la fermeture du 
puits, la courbe experimental est trade en echelle 
logarithmique en portant en abscisse les intervalles de temps A t 
et en ordonnee : 

20 t p + At 

At . AP« (15 > 

P 

t p represented la duree pendant laquelle la formation a Ste mise 
en production. V analyse des essais de puits pent alors etre 
faite par la comparison de cette courbe experimentale avec les 

25 courbes-type du graphe de la figure 2. 

La representation des courbes-type, avec en ordonnee P'jj-VS et 
en abscisse t^, est utilisable non seulement ? our D les 
formations souterraines homogenes mais egalement pour les 
formations non homogenes presentant par exemple une double 

30 P orosite. La figure 4 montre un exemple d' application a une 
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formation ay ant une double porosite. Dans ce cas ,le fluide 
produit par la formation est eontsnu daas la matrice, 
c'est-a-dire dans la roche composant la formation, et dans les 
interstices ou fissures contenues dans la matrice. On a done un 
5 systeme dans lequel le fluide. contenu dans la matrice s'ecoule 
d'abord dans les fissures avant de passer dans le puits. Le 
coefficient 03 caracterise le rapport du volume de fluide produit 
par les fissures au volume de fluide produit par le systeme total 
(matrice + fissure). Le coefficient X caracterise le retard de 

10 la matrice a produire le fluide dans les fissures par rapport a 
la production des fissures elles-memes. Le graphe de la figure A 
correspond a un modele theorique de formation ayant une double 
porosite. Sur ce graphe, on a represents, en traits pleins, les 
courbes-type correspondent au modele homogene, identiques a 

15 celles de la figure 2, en pointilles des courbes-type en 
choisissant comme index 

* C D 
1 -03 

et en semi-pointille des courbes-type en choisissant comme index 



03( 1 -03) 

Les courbes en pointilles represented l 1 equation : 

^^['■"'(■^ - 

25 Les courbes en semi-pointilles represented l 1 equation : 

• On a egalement represents par des points une courbe experimental 
typique caracterisant une formation a double porosite. 

30 L f utilisation du graphe de la figure A permet de determiner les 
valeurs des coefficients 0) et X , en plus des valeurs de 
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kh, C et S. On remarque que les courbes caracterisant le 
comportement d'un modele heterogene ont une forme tres marque* - P 
appliquant la methode selon la presente invention. 
La presente invention permet egalement de tracer sur un meme 
graphe theorique les courbes-type de la figure 2, P' t /C en 
fonction de t D /C D mais egalement les courbes-type P D en f one t ion 
de t D /C D decrites dans le brevet des Etats-Unis d'Amerique 
n" 4,328,705. La juxtaposition de ces deux series de courbes-type 
sur un meme graphe est montree sur la figure 5. On pent en effet 
realiser cette superposition sur un meme graphe car pour passer 
de P' D .t D /C D aux donnees experimentales qui sont A P'. At, il 
faut multiplier ces dernieres par un coefficient qui 1st donne 
par 1'equation (11). Pour passer de P D aux donnees experimentales 
_ AP, dans le cas des courbes-type du brevet cite, il f aut 
multiplier ces dernieres par le meme coefficient que 
precedemment. On pent done superposer les deux series de 
courbes-type et porter en ordonnee, avec une meme echelle, P et 
P VW Pour "tiliser le graphe theorique de la figure on 
utilise alors un meme graphe experimental comportant deux courbes 
representee en ordonnee les variations de pressionAP pour 1'une 
et AP VV C d P ° ur i'autre, At etant porte en abscisse pour les 
deux courbes. Le graphe combine de la figure 5 permet une 
comparison plus precise des deux courbes experimentales avec les 
courbes-type. 

25 La methode de determination des caracteristiques d'une formation 
souterraine qui vient d'etre decrite presente de nombreux 
avantages. Ainsi, 1'analyse des essais de puits pent etre 
effectuee a 1'aide d'un seul graphe, alors que les methodes de 
1'art anterieur font appel a un graphe general en echelle 

30 logarithmique en utilisant toutes les donnees experimentales et a 
un graphe specialise en echelle semi-logarithmique et ne prenant 
en compte qu'une partie des donnees experimentales. Du fait du 
comportement des modeles de systemes forma tion-puits" au debut et 
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en fin d'un essai de puits (temps courts et temps longs sur les 
graphes) qui se traduit par des lignes droites de pentes bien 
definies pour' les deux extremites "des courbes-type .le calage 
de la courbe experimental avec les courbes-type peut se faire 
5 sans ambiguite. La combinaison des courbes-type de l'art 
anterieur avec les courbes-type de la presente invention sur un 
meme graphe presente un avantage certain. De plus, la definition 
d'un nouveau temps, donne par l'equation (15), permet d' analyser 
les essais de puits effectues par la fermeture du puits. 

•0 II va sans dire que la presente invention ne se limite pas aux 
seuls modes de realisation qui ont ete decrits a titre d'exemples 
explicatifs mais nullement limitatifs. Ainsi, 1' evolution des 
valeurs de pression ou de la deriv€e des valeurs mesurees de 
pression peut etre comparee a Involution thiorique, calculee a 

«5 partir d'un modele de reservoir theorique, a l'aide de moyens 
informatiques tel qu'un ordinateur. 
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REVENDICATIONS 

1. Methode de determination des caracteristiques physiques d'un 
systeme forme d'un puits et d'une formation souterraine contenant 
un fluide et communiquant avec ledit puits, ladite formation 
presentant de l'effet parietal et/ou le fluide se comprimant ou 

5 se decomprimant dans le puits et ladite formation etant homogene 
ou heterogene, methode selon laquelle on provoque un changement 
de debit dudit fluide et on mesure une grandeur caracteristique 
de la pression P du fluide a des intervalles de temps successifs 
At, ladite methode etant caracterisee en ce que l'on compare, 
10 - d'une part, a partir d'un modele theorique de systeme de 
puits-reservoir, 1* evolution theorique du logarithme de la 
derivee P' D de la pression sans dimension en fonction du 
logarithme de t^, ladite derivee P' D etant par rapport a 
% ,C B' C D re P rgse ntant le temps sans dimension et le 
15 coefficient sans dimension de l'effet de compression ou de 
decompression du fluide dans le puits, avec 
- d' autre part, Involution experimental du logarithme de la 
derivee AP' de la pression en fonction du logarithme des 
intervalles de temps correspondants At, ladite derivee Ap' 
20 etant par rapport au temps t, 

et en ce que l'on determine de la comparison desdites evolutions 
theorique et experimental au moins une caracteristique du 
systeme puits-formation, choisie parmi le produit kh de la 
permeabilite k par l'epaisseur de ladite formation h, le 
25 coefficient C , 

2. Methode selon la revendication 1 , caracterisee en ce que la 
d'ite evolution theorique est celle du logarithme du produit 
P V t D /C D en fonction du logarithme de tp/Cp et ladite evolution 

30 experimental est celle du logarithme du produit AP'. At en 
fonction du logarithme de A t. 
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3. Methode selon la revendication 1 ou 2, caracterisee en ce 
que ladite evolution calculee est igalement reaction d'une 
grandeur caracteristique du produit C D e 2S et en ce qVon 
determine la valeur du coefficient S de l'effet parietal. 

4. Methode selon 1'une des revendications precedences 
caracterisee en ce que, lorsque la formation etudiee presente une 
double Porojit|, ladite^ Evolution theorique est en plus fonction 
des index | _ u e t - dans lesquelles X caracterise le 
retard de production de fluide par la roche de la formation 
souterraine compare a la production de fluide par les fissures de 
la formation souterraine et 0) repr€sente le rapport du volume de 
fluide produit par les dites fissures par le volume de fluide 
produit par le systeme total et en ce que l'on determine de la 
comparaison des evolutions theorique et experimental les valeurs 

15 de X et w. 

5. Methode selon l'une des revendications 2 a 4 caracterisee en 
ce que lorsque ledit changement de debit du fluide correspond a 
la fermeture du puits, on compare ladite evolution calculee avec 

20 Involution experimental du logarithme de 1' expression : 

t + A t 

-2 .At .AP« 

t 
P 

25 en fonction du logarithme des intervalles de temps A t, t etant 
le temps pendant lesquels le puits a ete mis en production! 

6. Methode selon les revendications 1 et 3 caracterisee en ce 
que l'on trace un graphe theorique de;.courbes r ty P e en coordonnees 

30 cartesiennes et en echelle logarithmique , ledit graphe 
representant ladite evolution theorique de la derivee P' , en 
fonction de tp/Cp. * D 
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7. Methode selon les revendications 4.et 6 caracterisees en ce 

qu'on trace en plus deux families dc c-i-S- -,- - 

aux index AV *D et D 

I -to co ( l -co) 

8. Methode selon l'une des revendications 6 et 7 caracterisee 
5 en ce qu'on trace une courbe experimental en coordonnees 

cartesiennes et avec la meme echelle logarithmique que ledit 
graphe theorique, ladite courbe experimental representant ladite 
evolution experimental de AP' en fonction de At, en ce qu'on 
fait correspondre ladite courbe experimental avec l'une des 
10 courbes-type dudit graphe theorique et en ce qu'on determine au 
moins l'une des caracteristiques kh, C, S, A et (o par le decalage 
des axes de coordonnees du graphe theorique et de la courbe 
experimental et par le choix de l'une des courbes-type. 

15 9. Methode selon la revendication 8, caracterisee en ce qu'on 
determine le coefficient C par le decalage des axes d'ordonnees 
de la courbe expSrimentale et du graphe theorique, kh par le 
decalage des axes d'abscisses de la courbe experimental et du 
graphe theorique et S, o et X par le choix de la courbe-type du 

20 graphe theorique correspondant a la courbe experimental . 

10. Methode selon les revendications 2 et 3 caracterise en ce 
qne 1'on trace un graphe theorique de courbes-type en coordonnes 
cartesiennes et en echelle logarithmique, ledit graphe 
25 representant ladite evolution theorique du produit P'.t /C , en 
fonction de t/C D D ' 



D' D" 



11. Methode selon les revendications 4 et 10 caracterisee en ce 
qu'on trace en plus deux families de courbes-type co.r^pondant 
30 aux index __D et AC D . 

1-to u(l-u) 
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12. Methode selon la revendication 10 ou 11 caracterisee en ce 
qu'on trace une courbe experimentale en coordonnees cartesiennes 
et avec la meme echelle logarithmique que ledit graphe* theorique , 
ladite courbe experimentale representant ladite evolution 

5 experimentale de produit AP'. At, en fonction de At, en ce 
qu'on fait correspondre ladite courbe experimentale avec l'une 
des courbes-type dudit graphe theorique et en ce qu'on determine 
au mo ins l f une des caracteristiques kh, C, S, A et W par le 
decalage des axes de coordonnees du graphe theorique et de la 

10 courbe experimentale et par le choix de l'une des courbes-type. 

13. Methode selon la revendication 12, caracterisee en ce qu f on 
determine le coefficient kh par le decalage des axes d'ordonnees 
de la courbe experimentale et du graphe theorique, C par le 
decalage des axes d 1 abscisses de la courbe experimentale et du 

15 graphe theorique et S par le choix de la courbe-type du graphe 
theorique correspondant a la courbe experimentale. 

1A. Graphes theoriques tels qu'obtenus par 1' application de la 
methode definie aux revendications 6 a 13. 
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